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INVESTIGATION OF THE EFFECT OF MATERIAL PROPERTIES ON 
COMPOSITE ABLATIVE MATERIAL BEHAVIOR 


by 

P. B. Cline and F. E. Schultz 


ABSTRACT 

The influence of material properties on the ablative performance of the silica phenolic, 
graphite phenolic and graphite epoxy materials was analytically investigated. Those pro- 
perties which have the greatest effect on the surface recession were established. The rea- 
sons for the property variation within a given material are stated along with the measured 
range of their values. 
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INTRODUCTION 


Prior to the development of analytical reaction kinetics ablation models, such as those 
developed under NASA contract NAS3-2566, the evaluation of ablative materials for use as 
rocket nozzle materials has been solely an empirical process. New or modified materials 
were made up and these materials were then test fired to determine whether they would 
perform better than previous materials. Although this approach does assess the gross ab- 
lative performance of a material in a given environment, it does not provide a detailed under- 
standing as to what particular parameters or properties of that material caused it to perform 
better or worse than some other materials. Recognizing these limitations of the empirical 
approach of evaluating rocket nozzle heat protection materials, several analytical methods 
have been derived for the theoretical prediction of the thermal performance of ablation ma- 
terials. One of these programs was the Reaction Kinetics Ablation Program (REKAP) which 
was developed by General Electric for NASA Lewis under contract NAS3-2566. This pro- 
gram is described in complete detail in Reference 20. 

The analytical approach of evaluating materials allows one to readily investigate the 
influence of the various material properties and environments on the thermal performance 
of the material. It was, therefore, the purpose of this contract, using the developed REKAP 
computer program, to establish the order of importance that the material and environment 
parameters have on the ablation performance (surface recession) of three materials. The 
materials considered within this study were silica phenolic, graphite phenolic and graphite 
epoxy. The steady state surface recession rates were calculated and the results are pre- 
sented as a function of the various parameters within this report. The first four material 
properties which had the greatest effect on the surface recession were then studied in detail. 
The property variations caused by processing, fabrication and design were investigated to 
establish a reasonable range for a comprehensive thermal investigation where the interplay 
of these four properties on the surface recession of each material was studied. 
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SUMMARY 


A three task analytical investigation was made to determine the influence of material 
property variations on the ablative performance of silica cloth/phenolic resin, graphite cloth/ 
phenolic resin and graphite cloth/epoxy resin materials as used for the thermal protection 
system of rocket nozzles. The thermal environment for which these materials were to be 
evaluated was that produced by the two propellant combinations, N„0 /Aerozine 50 and 
OF g/B 2^6* chamber pressure of 100 psia and operating at the combustion tem- 

perature corresponding to a 98% C* efficiency. The calculated surface recession at the 
nozzle throat was the ablative performance parameter of interest. The initial throat dia- 
meters were 1.5 inches and 7.82 inches. Only the environment resulting from NO/ 

Aerozine 50 was considered for the 7.82-inch throat nozzle. 2 4 

The thermal response of the materials, including the surface recession, was calculated 
using a Reaction Kinetics Ablation Computer Program. This program calculates the material 
temperature response including the thermochemical decomposition of the resin and the melting 
or oxidation of the reinforcing fibers or cloth. The ablation performance of a material de- 
pends on a number of material and environment properties. These properties are: 

Material Properties 


1. Char Density 

2. Virgin Material Density 

3. Char Thermal Conductivity 

4. Virgin Material Thermal Conductivity 

5. Char Specific Heat 

6. Virgin Material Specific Heat 

7 . Heat of Gasification 

8. Collision Frequency 

9. Activation Energy 

10. Heat of Melting or Heat of Reaction of Reinforcing Fibers 

11. Wall Emissivity 

12. Specific Heat of Ablation Gases 

13. Melting Temperature of Reinforcing Fibers 
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Environment Properties 


1. Heat Transfer Coefficient 

2. Recovery Temperature 

3. Nozzle Geometry 

The first phase of the study was to determine the amount of influence that each of the 
above properties had on the surface recession of the three materials. The material pro- 
perties presented in descending order of their relative effect on the ablative performance 
of each material are tabulated below: 



Silica Cloth/ 
Phenolic Resin 


Graphite Cloth/ 
Phenolic Resin 


Graphite Cloth/ 
Epoxy Resin 


1 . 

Melting Temperature 

.0288 

Surface Reaction Rate 
Constants 

.0064 

Surface Reaction 
Rate Constants 

.0082 

2. 

Virgin Plastic 
Density 

.007 

Char Density 

.002 

Char Density 

.0028 

3. 

Specific Heat (Char 
and Virgin Mat’l) 

.005 

Thermal Conductivity 
(Char and Virgin 
Material) 

.0017 

Thermal Conduct- 
ivity (Char and 
Virgin Material) 

.0018 

4. 

Activation Energy 

.003 

Virgin Plastic Density 

.00085 

Virgin Plastic 
Density 

.0011 

5. 

Char Density 

.002 

Collision Frequency 

.00065 

Activation Energy 

.0006 

6. 

Collision Frequency 

.001 

Specific Heat 
(Char and Virgin 
Material) 

.0004 

Specific Heat 
(Char and Virgin 
Material) 

.0004 

7. 

Heat of Gasification 

.0009 

Activation Energy 

.0002 

Specific Heat of 
Ablation Gases 

.0004 

. 8 ‘ 

Heat of Melting of 
Reinforcing Fibers 

.0008 

Specific Heat of 
Ablation Gases 

.0002 

Collision Fre- 
quency 

.0001 

9. 

Specific Heat of 
Ablation Gases 

.0008 

Wall Emissivity 

.0001 

Wall Emissivity 

.0001 

10. 

Thermal Conductivity 

.0008 

Heat of Gasification 

.0001 

Heat of Gasif- 

.0001 


(Char and Virgin Mat'l) ication 

11. Wall Emissivity .0005 
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The number alongside each property is the total variation in the surface recession rate 
(inches per second) over the range of interest of each parameter. 

Selecting the first four properties for each material, a detailed investigation of their 
variation, the reason for the variation, and the thermal response of the materials allowing 
each property to vary over its range, was made during Tasks n and HI of the program. The 
conclusions of the materials properly determination task (Task II) were that the material 
conductivities are more dependent on the lamination angle than on the resin or fiber content 
or the material density. The parameters found to have the major effect on the surface re- 
cession rate, such as the melting temperature of the silica fibers and the surface reaction 
rate constants of the carbonaceous char of the graphite cloth/phenolic resin or the graphite 
cloth/epoxy resin materials, are independent of the char or virgin plastic density, thermal 
conductivity, or activation energy. The activation energy is only a function of the resin 
material and not dependent on any of the other first three parameters. Therefore, it was 
concluded that each of the first four material properties which most affect the surface re- 
cession rate of these rocket nozzle materials are essentially independent of each other. 

Thus, during Task III the thermal performance (surface recession) of the three materials 
were analytically evaluated allowing each property to vary over its entire range regardless 
of the values of the other properties. 

The surface recession of the silica cloth/phenolic resin material was most influenced 
by the melting temperature of the silica fiber. The melting temperature range for the 
silica fiber extended from 3000°R to 4000°R. Increasing the melting temperature from 
3000°R to 4000°R decreased the recession rate by a factor of three. The virgin plastic 
density, specific heat and activation energy are the next most important material properties 
affecting the surface recession of the silica cloth/phenolic resin material. The combination 
of these properties in going from the minimum (88 lb/ft^) to the maximum (118 lb/ft^) value 
of virgin plastic density, from the maximum (75600 BTU/lb) to the minimum (21600 BTU/lb) 
value of activation energy and from the minimum (. 26 BTU/lb°R @530°R) to the maximum 
(.36 BTU/lb°R @ 530°R) value of specific heat increases the surface recession rate by 91 
percent. Therefore, it is the melting temperature of the silica fiber which has the greatest 
effect on the surface recession of the silica phenolic material. 

The surface reaction constants most influence the surface recession of both the graphite 
cloth/phenolic resin and the graphite cloth/epoxy resin materials. The surface reaction 
constants are defined by the following equation: 


q c 

m o K + K (h - h ) 

1 2 ' r w 

where mo is the mass loss from the surface of the material, is the convective heat flux 
received at the wall, Kj and K are the reaction constants, h^ is the recovery enthalpy; and 
h w is the enthalpy of the boundary layer gases evaluated at the wall temperature. The actual 
values of Ki and K 2 depend on chemical species in the boundary iayer in addition to the 
chemical composition of the surface of the material. The values of Kj and K 2 can be deter- 
mined either empirically, through the correlation of experimental data obtained from the 
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exposure of the material to the exhaust products of a rocket engine or, theoretically, through 
the use of a chemically reacting boundary layer program. Since the operating pressures of 
most rocket engines are relatively high ( > 10 atm) and the temperature of the exhaust pro- 
ducts is less than 8000°R, the graphite ablation is controlled by species diffusion within 
the boundary layer. Therefore, the surface mass loss may be calculated using a multi- 
component chemical model as developed by Scala and Gilbert (Reference 1). The results 
of these calculations are then correlated, yielding Kj and K 2 . For graphitic materials ex- 
posed to a turbulent air boundary layer, the constants Ki and K 2 are 4240 and 5.77 respec- 
tively. For the extremes of the chemical composition of the propellents considered in this 
study, the range of values for the constants were 1,000 to 12, 000 for Ki and 2 to 10 for K 2 . 

The surface recession rate of the graphite cloth/phenolic resin material increases from 
approximately 0.003 inches per second to approximately 0.024 to 0.027 inches per second 
in going from the maximum values of the surface reaction constants to their mi nimum value. 
Similarly for the graphite cloth/epoxy resin material, the surface recession rate increased 
from approximately 0.003 inches per second to 0.030 inches per second. The other three 
properties, char density, thermal conductivity and virgin plastic density, which most in- 
fluence the surface recession rate, have a considerably lesser effect on the surface recession. 
The combination of these properties in going from the maximum (82 lb/ft^) to the minimum 
(70 lb/ft^) value of char density, from the minimum (1.3 x 10"^ BTU/ft-sec°R @ 530 °R) to 
the maximum (11.3 x 10~^ BTU/ft-sec°R @530°R) value of thermal conductivity and from 
the maximum (95 lb/ft^) to the minimum (85 lb/ft^) value of virgin plastic density causes the 
surface recession rate of the graphite cloth/phenolic resin material to increase only 44 per 
cent. Similarly the surface recession rate for graphite cloth/epoxy resin material increases 
only 39 per cent in allowing the char density to decrease from 76 to 64 lb/ft^, the thermal 
conductivity to increase from its minimum to its maximum value (1. 2 x 10“^ to 11.4 x 10 - ^ 
BTU/ft-sec°R @ 530°R) and the virgin plastic density to decrease from 95 lb/ft^ to 85 lb/ft**. 
Therefore, to minimize the surface recession of graphitic materials, additives or methods 
of protecting the nozzle walls from the reactive chemical species in the boundary are desirable. 
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DISCUSSION 


The analytical investigation of the influence of the material properties and the external 
environment on the ablative performance of silica phenolic, graphite phenolic and graphite 
epoxy was divided into three phases or tasks. The first task was a screening investigation 
to determine the relative importance of each parameter in the ablation or surface recession 
process. The second task was to determine the range of variation of the four most impor- 
tant material properties affecting the surface recession. The third task was a detailed 
study of the ablation performance as these four parameters were allowed to change. The 
analytical evaluation of the material performance was accomplished by using the Reaction 
Kinetics Ablation Program (REKAP) developed on NASA Lewis Contract NAS3-2566. A 
mathematical description of the program is given in Appendix A. 

The screening investigation involved evaluating the steady surface recession rate of 
the three materials for each of several variables. The variables included material pro- 
perties and environments. The variables considered were: 

Material Properties 

1. Char Density 

2. Virgin Plastic Density 

3. Char Thermal Conductivity 

4. Virgin Material Thermal Conductivity 

5. Char Specific Heat 

6. Virgin Material Specific Heat 

7. Heat of Gasification 

8. Collision Frequency 

9. Activation Energy 

10. Heat of Melting or Heat of Reaction of Reinforcing Fibers 

11. Wall Emissivity 

12. Specific Heat of Ablation Gases 

13. Melting Temperature of Reinforcing Fibers 
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Environment Parameters 


1. Heat Transfer Coefficient 

2. Recovery Temperature 

3. Nozzle Geometry 

The ranges of the material properties used for the screening process are tabulated in 
Table 1. The nominal material property values are those which are commonly accepted as 
being the average properties for these classes of materials. The range of each variable 
is based on numerous ground test results and the judgement of material manufacturers of 
what future modification in the material formulation would do to the material properties. 

The environment parameters were those calculated for the nozzle throat, assuming the pro- 
pellant combinations to be either N204/Aerozine 50 or OF2/B2H6, throat diameters of 1.2 
and 7. 82 inches, a chamber pressure of 100 psia and a combustion temperature (or re- 
covery temperature) corresponding to a C* efficiency of 96 per cent. Maximum environ- 
ment parameters for the graphite phenolic and graphite epoxy materials resulted from the 
OF2/B2H6 propellant combination but are limited to the 1. 2 inch diameter throat. The 
range of environments for the silica phenolic material results from considering both the 
1.2 and 7. 82 inch diameter throat nozzle using the ^O^/Aerozine 50 propellant combina- 
tion only. The procedure followed during the screening effort was to vary each of the para- 
meters given in Table 1 individually while holding the rest of the parameters at their 
nominal values. Tabulated in Table 2 in descending order are the parameters and their 
relative importance in effecting the surface recession of the three materials. The relative 
importance parameter is the change in surface recession rate as the parameter is varied 
over its range. The actual surface recessions as functions of the various parameters, are 
presented in Figures 6 through 44. The thermal conductivities and specific heats are pre- 
sented as ratios (the value under investigation divided by the nominal value) rather than the 
absolute value, since in this way the temperature variation was eliminated. 


The environment parameters (film coefficient and recovery temperature) appeared 
among the first four parameters which most affect the surface recession rate of all mater- 
ials. Since these parameters were independent of material, they were not investigated any 
further. The nominal environment conditions corresponding to the propellants and nozzle 
sizes were used during phase three of the investigation. The four material properties 
selected for the detailed investigation were: 


Silica Cloth/ 
Phenolic Resin 

Melting Temperature of 
Reinforcing Fibers 

Virgin Plastic Density 

Specific Heat (Solid) 

Activation Energy 


Graphite Cloth/ 
Phenolic Resin 

Surface Reaction Constants 
Kj and K2 

Char Density 

Thermal Conductivity 

Virgin Plastic Density 


Graphite Cloth/ 
Epoxy Resin 

Surface Reaction 
Constants Ki and K2 

Char Density 

Thermal Conductivity 

Virgin Plastic Density 


8 


The surface reaction constants Kj and K 2 are dependent on the composition of the 
boundary layer gases and on the chemical composition of the ablative material. Since these 
constants describe a single surface reaction rate, they were grouped together and were 
considered to be a single value during Phase HI. 


Tabulated in Table 3 in descending order are the material properties which most af- 
fect the steady-state mass loss rates. The "Relative Importance Parameter" columns 
within this table give the mass loss rate differences between that occuring at the maximum 
and minimum values of each of the properties. Mass loss rates are shown in Figures 45 
through 80. 

Tabulated in Table 4 is the amount of internal degradation, or increase in char thick- 
ness, occurring during the nozzle cool-down period. The nominal degradation during noz- 
zle cool-down of the Silica Cloth/Phenolic Resin material is approximately 1/6 that of the 
other two materials. This is due to the fact that less thermal energy is stored in the char 
material. Similarly, it can be seen that when the thermal or structural properties are such 
that very little char is formed (for example, the low melting temperature of the reinforcing 
fiber condition), the internal degradation during cool-down is minimal. 

The Phase II portion was a study to characterize the variation in thermal properties 
resulting from processing and fabrication techniques for silica-phenolic, graphite- 
phenolic, and graphite-epoxy materials. This effort was primarily accomplished by con- 
ducting an extensive literature search of existing data generated by both government and 
industry sponsored efforts. Certain prominent trends and peculiarities were discernible 
m the variations m thermal conductivity and specific heat coefficients due to change in 
lamination angle, resin content, density, and lot-to-lot variation. 

Table 5 identifies the graphite -phenolic thermal conductivity data plotted in Figure 
82 with respect to lamination angle, resin content, density, and source of information. 
Figure 82 illustrates, in general, that an increase in the resin content of the graphite- 
phenolic laminates, keeping lamination angle constant, decreases the thermal conductivity. 
For example, an increase in resin from 30 percent to 50 percent decreases both the across 
lamina and with lamina direction (Figure 81) thermal conductivity coefficients for a partic- 
ular resin-cloth composite by approximately 43 percent (Figure 82, curves 13 and 16) and 
25 percent (Figure 82, curves 1 and 6) respectively (See Reference 13). Also, a decrease 
in thermal conductivity with decreasing density was observed for the across lamina speci- 
mens. This trend is anticipated since an increase in resin content would be expected to re- 
su t in a graphite -phenolic composite of decreased density and lower interlamina thermal 
conductance. This dependence of thermal conductivity upon changes in resin content and/ 
or density may well explain the variation between lots of graphite-phenolic. 

It was observed that the width of the scatter bands for the thermal conductivity data of 
the across lamina and with lamina specimens was comparable. However, the values of 
thermal conductivity for the with lamina specimens are considerably higher than those for 
the across lamina. This is true because a more direct path for heat flow exists in the warp 
or fill direction compared to the across lamina path with its discontinuous resin-cloth 
layered configuration. Also, for the "with lamina direction" one observes higher thermal 
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conductivity coefficients for the warp direction whose fibers are more direct and usually 
more numerous than for the woven fill direction which normally contains less fibers per 
unit length. For example, a 23 percent difference in thermal conductivity coefficients be- 
tween the warp and fill directions was observed for two graphite -phenolic specimens of 
identical composition (See Figure 82). This points out the hazard in simply defining a single 
thermal conductivity function for this type of material. 

Figure 84 shows a slight variation in specific heat for the graphite -phenolic materials. 
This is due to the fact that specific heat does not appear to be as strong a function as ther- 
mal conductivity of resin content, density, or lot-to-lot variation. It can also be seen that, 
as expected, lamination angle has no effect upon specific heat. 

Based upon the data listed in Tables 5 and 6 and plotted in Figures 82 and 83, nominal 
curves for the thermal conductivity and specific heat coefficients of graphite -phenolic are 
shown with a tolerance band in Figures 84 and 85, respectively. Due to the scarcity of char 
data, the tolerances placed on the virgin material thermal conductivity and specific heat 
nominal curves were extended to the char portion of the curves. The tolerance bands for 
graphite -phenolic composites, as well as the other subject materials, are wider than those 
which appeared in the preliminary release as a result of the evaluation of additional data. 

The previous trends observed for graphite-phenolic are not discernible for silica- 
phenolic. In fact, the trend is reversed for one composite grouping where a lower density 
is recorded in both the across and with lamina directions for the lower resin content 
composites (See Tables 7 and 8). This suggests that the creation of voids may have been 
associated with the reduction in resin content. The presence of high void content would 
be expected to yield lower thermal conductivity coefficients. 

Referring to Figure 86, it can be seen that the with lamina materials exhibit higher 
thermal conductivity coefficients than the across lamina materials. Also, the spread in 
thermal conductivity data for the with lamina composites is less than that for the across 
lamina composites. 

A slight variation in the specific heat of the silica-phenolic composites can be seen in 
Figure 87. The specific heat coefficients of these composites are lower than those of 
graphite -phenolic and graphite -epoxy. Since the silica cloth has a lower specific heat than 
graphite cloth, it is not surprising that this trend is retained when each system is intro- 
duced as a woven cloth in a resin binder. It should also be noted that the specific heat of 
the phenolic resin is significantly higher than that of the reinforcements of Silica and graph- 
ite resulting in composite specific heat coefficients lying between those of the resin and the 
reinforcement. 

Nominal curves for the thermal conductivity and specific heat coefficients of silica- 
phenolic are shown with a tolerance band in Figures 88 and 89, and are based on the data 
presented in Figures 86 and 87. It should be noted that the tolerances placed on the virgin 
material thermal conductivity and specific heat curves were extended to the char portion 
of the curve because of the lack of char data. In fact, the thermal conductivity and specific 
heat coefficients of fused amorphous silica at elevated temperatures served as the basis for 
the high temperature trend (See Figures 86 and 87). 
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For the graphite -epoxy system, experimental thermal properties data is very limited. 
Therefore, the nominal curves of thermal conductivity and specific heat coefficients of 
graphite -epoxy presented (in Figures 90 and 85, respectively) are based on the trends 
established for graphite -phenolic because of the similarity that exists between the thermal 
conductivity and specific heat coefficients of the phenolic and epoxy resin systems, with one 
exception (Reference 10). This exception occurs in the region of resin decomposition. In 
this region the thermal conductivity of graphite -epoxy is shown with a steeper slope than 
that of graphite -phenolic. This arises because epoxy resins lose approximately 80 percent 
of their initial weight in this region of decomposition, whereas the phenolic resins lose 
about 47 percent of their initial weight in this region which extends over a wider tempera- 
ture range. In addition it should be noted that the epoxy systems have poor thermal stab- 
ility at high temperatures, and therefore, are not desirable for rocket nozzle design. 

A fiber reinforced plastic, such as discussed in the context of this report, is normally 
fabricated from prepreg materials. This prepreg is a web material such as carbon cloth, 
refrasil cloth, etc. , impregnated with a resin. The resin is advanced to a dry stage by 
partial curing. This operation is continuous, and is accomplished by passing the web through 
a pan of resin, then metering through squeeze rolls, followed by resin advancement into 
temperature controlled ovens or towers. From these prepreg materials, various shapes 
are laid up in layer form utilizing predesigned patterns. These materials are then fully 
cured in molds under heat and pressure to attain the desired configuration. In most cases, 
additional machining techniques are required to complete the fabrication procedure. 

The pressure employed in molding is dependent upon the resin system employed. Cure 
of epoxy type resins is by addition; i. e. no volatile products are produced or evolved in the 
cure mechanism. Thus relatively low pressures are required in laminating. The phenolic 
resins cure by condensation; i. e. , water is a by-product of the cure mechanism. Thus to 
maintain the desired high density of the molded part, high pressure, in excess of 1000 psi, 
is required in molding these materials. 

The curing temperature is dependent upon the resin system. Temperatures normally 
applied are in the 300°F range. Normal practice is to step-wise attain the maximum tem- 
perature and to hold sufficiently long enough to insure thorough resin cure. This time factor 
also is dependent upon the resin system employed. 

Depending on several factors such as shape, size, materials, etc. , the part may be 
fabricated between platens, by hydraulic press or in a hydroclave, (a hydroclave is a 
pressure container filled with water; the sealed, immersed laminate is cured by the hot 
water under hydraulic pressure). 

The amount of resin, the degree of resin advancement, and the volatile content of the 
prepreg material is carefully controlled prior to laminating. A resin advanced too far will 
not flow and knit in the laminating operation. A resin with excess flow will squeeze out, 
and result in a resin starved laminate. With phenolic resins, a critical level of volatile is 
essential to plasticize the resin, and obtain the desired resin flow. Excessive volatiles 
distort particular laminate properties - dielectric, ablation, etc. 
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When a load is applied to a fiber reinforced plastic material, a complex stress pattern 
is produced within the material depending upon the orientation and proportion of the fiber- 
matrix as well as strength and adhesive properties of the resinous component. In these 
composite plastic structures, the resin supports the reinforcement under hydrostatic pres- 
sure, and the reinforcement continues to function as a load-carrying member up to the 
point at which this pressure decreases as the resin yields, following which the composite 
fails. A common mode of tensile failure is a delamination between plies and across plies. 

The limiting strength of current unreinforced resin systems appears to be about 15,000 
psi in tension and 25, 000 psi in compression. Thus, to attain improved physical properties, 
reinforcement materials are employed with the resinous adhesives. The choice of these 
materials are usually selected by trade-off studies, and are strictly dependent upon the de- 
mands imposed by the end application. 

From a consideration of the variables involved in making reinforced shapes, it is 
apparent that a process history and web orientation be established before acceptable physical 
properties can be determined in a uniform and realistic manner. Some variation is to be 
expected when comparing properties provided by several techniques and on parts prepared 
using non-uniform processing. 

The manufacture of sophisticated items such as an ablative rocket nozzle requires the 
utlization of all of the art and science at the disposal of the laminating industry. Due to the 
non-homogeneity of the product, the resulting properties are directionally dependent and 
subject to variations caused by the processing techniques. A specific resin fiber system 
is seldom processed in the same manner by individual vendors which in turn potentially 
introduces further property variations within the laminates. The variations in vendor pro- 
cessing are, most often, due to capability of fabricating equipment and manufacturing philos- 
ophy. In general, the methods employed are company proprietary and are net for general 
publication. 

The properties of laminates will vary with different levels of resin content. Of course, 
many additional factors come into play. The variation in resin content, or density of the 
laminate, however, is a very important processing factor and much effort is extended to 
control the ratios of reinforcement to resin within the laminate, hi addition to the propor- 
tions of these two ingredients, there are other considerations such as the relative location 
of the resin on and in the web, as well as the degree of resin advancement prior to molding. 
Additional factors which influence final properties at this point of fabrication include age of 
material prior to laminating, residual volatiles, etc. 

The mechanism of impregnating and wetting the fibers in the reinforcing bundle is to 
utilize a solvent solution of the resin. The solvents employed with the systems discussed 
herein are normally low boiling materials such as alcohol, toluene, etc. The solution may 
contain additional proprietary ingredients such as wetting agents, surface activating agents, 
fillers, catalysts, etc. 
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The surfaces of filamentary and fibrous materials are not always receptive to chem- 
ical bonding to matrices, so that untreated fiber surfaces are rarely used in the prepara- 
tion of reinforced composites. A specific example may be found in the use of finishes of 
various types of surfaces of glass fibers. Silane or chrome finishes with appropriate or- 
ganic groups in the molecule are widely used to form a ,r bridge'' between the glass surface 
and the resin; the silane portion of the bridging agent reacts with the O-Si-O molecules in the 
glass, and the organic portion is free to react with the resinous binder. In some instances, 
this bridging agent is incorporated in the impregnating solution and in others, it is applied 
by a separate operation to the glass surface prior to resin impregnation. The use of this 
agent is not reflected in the density or resin content of the final laminate; it does, however, 
display its absence by changes in mechanical properties upon aging. 

In some instances, such as with carbon, improved fiber wetting by the phenolic resin 
is realized with traces of proprietary wetting agents added to the resin solution. This item 
promotes greater uniformity in resin distribution and thus enhances the uniformity of the 
laminated product. 

A factor in considering carbon webs initially was the wide variation of the active sur- 
face area within a given batch of carbon cloth. This variation resulted in poor impregna- 
tion and excessive localized deposits of resin on the web surface. This condition very 
likely distorted physical properties. The vendors of carbon cloth today appear to have this 
problem under satisfactory control. 

Test data developed in the evaluation of carbon phenolic composites in rocket engines 
indicated the need of a carbon filler in addition to the carbon web within the fabricated 
parts. These data indicated improved ablation resistance, thus, the same materials are 
durrently employed in carbon base heat shields. Studies are under way to establish the 
need of these fillers. The presence of these carbon fillers within the resin solution re- 
duces or retards resin impregnation into the fiber bundle, deposits the particles onto the 
web surface and results in a general lowering of the potential available interply adhesive 
bond strength, which thus results in lower flexural strength of the finished laminates. 

Of the systems considered in this evaluation, only the epoxy formulations incorporate 
a catalyst component. The phenolic resins cure by condensation reactions under heat and 
high pressure; the epoxy resins cure by addition in the presence of a catalytic agent under 
heat and relatively moderate pressures. In general, epoxy systems may have a greater 
density spread due to the lower fabricating pressures. 

The observed density variations of 35 phenolic carbon and 48 phenolic refrasil items 
recently purchased to a given density specification of 92.5 lb/ft 3 for phenolic carbon and 
105 lb/ft 3 for phenolic refrasil are as follows: 
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Density of Laminates (lb/ ft^) 


Phenolic Carbon Phenolic Refrasil 


Min. 90. 5 99. 5 

Max. 93 110 

Average 92.5 105 

Range -2, + 1/2 ± 5 

O 

The refrasil fibers have a density of 131 lb/ft and the carbon 106 lb/ft 3 . Data on the 
graphite phenolic was not available, but it is expected that its density variation would be 
similar to that of the carbon base laminates. Data on epoxy composite formulations was 
not obtained. 

The narrow spread in density range for the carbon phenolic system is indicative of the 
excellent control of the processing variables. The spread in density range for phenolic 
refrasil, though broader, is reasonable. 


The relationship of density of a laminate to some physical properties has been estab- 
lished. Graphic representation of these relationships, assuming proper fiber treatment, 
displays optimum or maximum values over a relatively narrow spread in resin content. 
Resin content by itself, however, is not indicative of the material performance. The dis- 
tribution of the resin within the prepreg and the production and fabrication procedures of 
the laminate combined, impart the final system properties. 

With a given resin system, differences exist in the impregnating characteristic of 
various webs. Thus the ability of the resin to wet through and impregnate the center of the 
fiber bundle contributes to the variability of the product density, etc. 

For example, the fibrous bundles, be they carbon, refrasil, asbestos, etc. , in each 
case require a certain amount of resin to restrain their relative motion by filling in the 
natural voids within the bundle. In turn, the voids between bundles, and likewise, the 
voids between plies, etc. , each require a specific resin content which, commensurate with 
optimum fabrication procedures, reflect maximum or ultimate performance. Thus a spe- 
cific resin content in a given laminate may yield a maximum compressive strength, whereas 
a different level of resin content may be required to realize the ultimate in interlaminar 
bond strength. Obviously, in these instances, a trade-off is essential to comply with the 
overall desired end properties of the laminate. 

Thus in rocket engines, the mechanical ablative performance is determined by the 
selection of raw materials, and each succeeding step in combining these materials and 
processing them into finished hardware. Control of the potential variables in overall 
processing of the materials is essential to realize ultimate quality parts. 
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Based on the results of the material property investigation (Phase II), the variations 
in the thermal conductivities were more dependent on lamination angles than on the resin 
or fiber content. The parameters which were shown to have a major effect on the surface 
recession rate, such as melting temperature of the silica fibers and the surface reaction 
constants of the carbonaceous char, were independent of the char or virgin plastic densi- 
ties, thermal conductivity, or activation energy. The activation energy was a function only 
of the resin material and not dependent on the other three parameters. Therefore, it was 
concluded that each of the properties is independent of the others. The material perfor- 
mance during the detailed investigation of Phase III was calculated allowing each of the 
properties to vary over their entire range. 

The ablative performance during Phase III was calculated using the information given 
in Tables 9, 10, 11 and 12 and in Figures 84, 85, 88, 89 and 90. The surface recession 
rates for silica cloth/phenolic resin are plotted in Figures 6 through 23 for the range of 
material and environment properties given in Table 1. Figures 91 through 99 show the 
recession rates obtained for the conditions existing within the 1.2 inch throat diameter 
nozzle, while Figures 100 through 108 give the recession rates for a 7.82 inch throat dia- 
meter nozzle. Figures 91 through 108 indicate the following descending order of property 
influence on the recession rate of silica cloth/phenolic resin A, for the property variations 
given in Table 1: silica fiber melting temperature, virgin plastic density, specific heat, 
activation energy. The specific heat is presented as the ratio of the value used to the nomi- 
nal value. The recession rates for the 1.2 inch throat diameter nozzle were slightly higher 
than those for the 7 . 82 inch throat diameter nozzle, which is expected since the convective 
heat transfer coefficient is higher for the smaller throat diameter (0.294 BTU/ft. 2 sec °F 
as compared to 0.223 BTU/ft. 2 sec °F for the larger throat diameter). 

The surface recession rates for the graphite cloth/phenolic resin and for the graphite 
cloth/epoxy resin materials are shown in Figures 109 through 117 and in Figures 118 
through 126 respectively. For both of these materials, the property which has the strongest 
influence on surface recession rate is the reaction rate (constants K-^ and K 2 ), followed (in 
descending order) by char density, thermal conductivity and virgin plastic density. The 
thermal conductivity is divided by the nominal value in order to eliminate the effect of the 
variation of conductivity with temperature (Figures 115-117 and 124-126). 
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TABLE 1. PROPERTY VARIATION FOR SCREENING INVESTIGATION (CONT) 


GRAPHITE CLOTH/EPOXY RESIN 

Max 

0.75 

8160 

» 

Nominal 

0.425 

4990 

7.82 

Min 

0.10 

3500 

1.5 

GRAPHITE CLOTH/PHENOLIC RESIN 

Max 

0.7 

8160 

os 

Nominal 

0.425 

4990 

7.82 

Min 

0,10 

3500 

1.5 

SILICA CLOTH/PHENOLIC RESIN 

Max 

0.47 

5630 

00 

Nominal 

0.294 

4565 

7.82 

Min 

0.118 

3500 

1.5 

PROPERTY 

Environment Parameters 

1. Heat Transfer 
Coefficient 
(Btu/ft^-seo^) 

2. Recovery 
Temperature (°R) 

3. Nozzle Geometry 
(Throat Diameter 
In Xnohes) 


17 











Ills 

® a § 5 

* * 


CO i— I 

to 

o o 
o o 


rH <X> 
rH O 
O O 
O O 


o o 
o o 
o o 


o o 
o o 
o o 


w v»n 

S S& 

8 ** 


3 r?> £> 


§ Q ff 

■5 o 2 


r * 

O cd 

-s S 


^ CQ 

ng (D 

S CQ 

O <5 

CQ 0 


cd aj 

CD CD 

K K 


£ Sd 

CO CO 

co cd 

*s o 


F*h CQ O O H > < 


o o 

CD 0 

a a 

CQ CQ 


O <D 

U £ S 


® u 
o 0 ) 

a 4-> 

Cd CD 

O 

in 


CM 

i— i 


t- 

in 

00 

m 

CO 


CM 

CM 

rH 

rH 

t 5 8 

© 

oo 

CO 


CO 

CM 

rH 

o 

o 

o 

o 

o 

o 

O 

M M 

o cd 

a u 

CQ 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d cd 

*9 ft 

d 

d 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

ft 

o 

ft 

o 

ft 

o 

d 


q 9 ^ £> 

5 J ^ W ^ 
«B <2 & co 
© Ph d 

» si a 

a 1 1 | 

•'— 1 3 O rd 

Ph CQ O O 


^3 CQ 

> d 

r* © 

I ° 

-© o 


O rt 

*3 ^ 
a q 


o 3 

Cl 

CD O 
d CQ 
CJ 1 v ^' 
© 

cd 

Pm ,2 


o o 

*t -4 ;r-( 

.2 sa 


& 

© 

CQ 

d 


c 

o 

•r-< 

2 

o 

& 
• 1— 1 

cd 

o 

d 

W 

td 

© 

CQ 

3 

CQ 

d 

o 

a 

CQ 

•pH 

ci 

cd 

O 

•pH 

td 

o 

a 

s 

w 

MH 

o 

> 

•pH 



•pH 

•M 

o 

< 

O 

© 

a 

CQ 

Id 

td 

© 

a 


5og< 

® §* c 3 .a 

Ph 3 cl w 


00 

LO 

00 

00 

m 

co 

CM 

CM 

CM 

O 

o 

o 

ft 


CO CM 

o o 
o o 


O 00 

i— I o o 

o o o 

o o o 


oo oo in 

o o o 

o o o 

o o o 


d co cd 

t* Sh U 

u J 8 8 , 

& e a 

s bD £ 

0 ) g h 

H -S >, 

CJ f, 

» *a ® 

•S .2 > 

a c © 

,2 ® o 

S 3 « « 


£ W 

W cd 


S >> 

O ho 

CQ ^ 

w 2 

+2 d 

Cd H 

/vv 


‘5 .s a 

2 o £ 

M <1 O 


o .2 w S 

•rd *“• 2 

cd © <5 

§ .2 a 2 . 

s o -S 

2 a bo 2 


cd a 

o > -g 

^ 4-1 

o o o 

*+5 *3 

cd cd -d 

a) © © 


vj co 

O CO 


o cd 

3 ft 

a> 2 


18 



ABLATIVE MATERIAL DEGRADATION DURING NOZZLE COOL-DOWN 
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TABLE 5. GRAPHITE CLOTH/PHENOLIC RESIN - THERMAL CONDUCTIVITY 
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Here, Wj is the net amount of specie i produced per unit volume per unit time. Note that 
W. includes the formation of the species from the unreacted solid as well as any further gas 
phase or gas-solid phase reactions that might occur. 


The surface integral involved in (13) can be transformed into a volume integral by means 
of the Divergence Theorem*: 


/ ^ aA= / ’'p,/v+v d \ 


dV 


The order of integration and differentiation can be interchanged* so that: 

i / p i dV= /ir dV 

V V 

Substituting these relations into (13) yields: 


/ ir + 7 -pi(* + v d .)-w, 

v L ' 1 


dV = 0 


Since the volume is arbitrary the integrand must be identically zero. Thus, the species 
continuity equation is: 


* + 7 -*i( 1 ?+ v) = w i 


(14) 


Summing this equation over all gaseous species and noting that 


£ "i v d. = 0: £ p i = V £ V w g 


• o 


results in the continuity equation for the gas: 

dp 

IT + v * P V* W 
9t g g 


(15) 


*It is assumed that all functions are continuous and continuously differentiable and that the 
region is simply connected (Reference 6). 
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Now, gas phase reactions do not change the total mass of gas present; rather, they redis- 
tribute the species. Therefore, W = ? W. is the total rate at which gas is being produced 

o 1 1 

by the decomposition of the unreacted material and by gas-solid phase reactions. 

The continuity equations for the solid species are: 


EL 

at 

= w 

p 

(16) 

ap c 

at 

= w 

c 

(17) 


W is the rate of depletion of the unreacted material due to decomposition. For charring 
ablation materials, the decomposition is irreversible and the rate at which it proceeds is 
generally limited by chemical kinetics. W p is deduced from TGA (thermogravimetric 
analysis) experiments and is often expressed analytically as a single nth order reaction 
with an Arrhenius "rate constant". 


W = - A (pj n 

p ' y P' 

—) 

RT J 

W c is the rate at which char is formed from the decomposition of the unreacted material 
(generally a known fraction of W p ) plus gas-solid phase reactions. There is no overall 
production of mass so that: 

W + w + W =0 ( 18 ) 

peg v ’ 

It is useful to separate out the mass production rates due to the decomposition, the gas 
phase reactions, and the gas-solid phase reactions. This can be done by introducing some 
new quantities. 

. . . »» 

W = - f w + w 
c c p c 

W = - (1 - f ) W + W 
g c P g 

Here f c denotes the fraction of unreacted material which forms char (not necessarily con- 
stant) and the superscript double prime denotes gas-solid phase reactions only. The first 
part of these equations state that char and gas are produced from the decomposition of the 
unreacted material, while the second part accounts for additional formation due to gas-solid 
phase reactions. Note that 


( 19 ) 

( 20 ) 


A = A exp 
o 
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( 21 ) 


• it , v 

W = - W 
g c 


which follows from (18). 

Finally, the species continuity equation can be expressed in terms of gas phase reactions 
only. Using p and the chain rule on (14) yields: 

o 


i at 


+ V * o v 


. +p g [v + v ' VK i + 7 • ( p g v d. K i ) = w, 


The first term in brackets equals W by virtue of (15) and so: 

o 

"3K. 1 / 

p g [^r + v ' VK iJ 


Where: 

W. = W. - K. W 
i i i g 


(23) 


Wj is the net rate of production of the ith species minus the amount of the ith species formed 
by the decomposition of the unreacted material plus gas-solid phase reactions. Consequently, 
w i is the net ra te of production due to gas phase reactions only. Note that: 



• f 

In general, the are functions of temperature, pressure and composition and are deter- 
mined from a knowledge of the exact chemical reactions (and these rates) which occur. For 
very slow reactions in the gas phase (i.e. , "frozen-flow"), the = 0. For very fast 
reactions, the flow will be in local thermochemical equilibrium and the Wj' are determined 
by imposing constraints on the composition (i.e., equilibrium "constants"). 

Energy Equation 


The energy equation is derived by applying the First Law of Thermodynamics to a station- 
ary control volume within the material. This means that the time rate of change of the total 
energy within the volume equals the rate at which energy is transported into the volume 
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minus the rate at which energy is being convected out plus the rate at which work is being 
done on the volume. The mathematical expression is 


_d 

dt 




Q • n dA 


V dA 


(24) 


The energy flux vector Q may be expressed in terms of contributions due to heat conduction, 
thermal radiation and diffusion (References 2-4): 

5 = VV^ p iV h i 
1 1 


It now remains to relate the heat flux vectors and the surface force per unit area to the 
variables of the problem. Since the solid and the gas are in intimate contact, it is assumed 
that the temperature of the gas equals that of the solid. The conduction heat flux vector is 
approximately linearly dependent upon the temperature gradients. For an isotropic material* 
this implies 

"q = - K V T 
c 

which is Fourier's Law. 

For an isotropic material, the heat flux depends upon temperature gradients through a 
second order conductivity tehsor. In rectangular cartesion coordinates, this is (p. 38, 
Reference 7): 



5T 

5x j 


e 

i 


(25) 


For example, the ablation material may be somewhat anisotropic due to fiber-type fillers 
in the solid material or because of the changes in composition. ** 


♦Isotropic material - medium whose structure and properties in the neighborhood of any 
point are the same relative to all directions through the point (p. 6, Ref. 7) 

**An excellent example of an anisotropic material of current interest is pyrolytic graphite 
but since it does not decompose in depth it is not pertinent to the present problem. 
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For the purposes of the present analysis, the material will be considered isotropic 
although it is noted that it would be easy to include (25) in the analysis should sufficient data 
be available to justify it. Thus, the conduction heat flux vector is related to the temperature 
by: 

■q = - K V T (26) 

c 

Note that the conductivity will be a weighed average of the conductivities of all species that 
are present, both solid and gas. 


In general, the radiation heat flux vector accounts for the net effect of emission, absorp- 
tion and scattering of thermal radiation of all wavelengths within the material. It is usually 
assumed that scattering is negligible, the material is isotropic and that the optical proper- 
ties do not depend on the wavelength. Even with these drastic assumptions, the calculation 
of qp is quite complex. Thus, for practical calculations, q R is usually neglected and the 
transport of thermal radiation is approximately accounted for by an increase in thermal 
conductivity with temperature. 


The surface force per unit area can be related to the stresses by considering the forces 
acting on a small tetrahedron (p. 101, Reference 8). 


-i. 

P = n 


T 


n is an outward unit normal and T is a second order stress tensor whose components are 



Note that this is a dot product of a vector with a tensor and the result is a vector which is 
different than n, both in magnitude and direction. Using indicial notation, the surface force 
would be expressed as: 


P = a. 

u 


n. e. 
i 1 


The required work term is: 


p . V = ( V • T ) • n 

The stress tensor may be separated into a hydrostatic pressure component (a scalar) 
and a viscous stress tensor. 


T = - P + r 


For a linear isotropic fluid, the viscous stresses are linearly related to the velocity 
gradients. For the purposes of this analysis, the viscous stresses will be retained in the 
general form of T . The final form of the work term is then: 

P • V s - (p V) • n + (V> T) • n 


A-ll 



Substituting (26) and (27) into (24) and following exactly the same procedure as with the 
species continuity equation results in the differential energy equation: 



The various terms of the energy equation can be expanded and rearranged to a more co 
convenient form. Adding the term + V • p V - to both sides of (28) and using Dalton's Law 
gives: 



and using the chain rule, (29) can be expanded to: 
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Combining these relations, the final form of the second term of (31) is: 

[• * •] = ' W ft 1 " U h + f e - e l + E W. ? h. - w" (h - e ) 

LJ PL cgcc vpj . i i c ' g v' 


( 32 ) 


The eighth term of (31) can be expanded to: 


[• * •] 


V • V 


dp 


g 


+ V . P V 
d t 'g 


+ P 


g 




Using (15), the final form of the eighth term is: 


V . V • 

[•••)■ -r- w , * », 


(V) 


Substituting (32) and (33) into (31) yields a final form of the energy equation. 
Storage 


(' 


pc + p C +pc \ — 
P v p M c v p p g pj at 


(33) 


- W 


Decomposition 


(l-f)h + f e -e 
c g c c p 


Cracking Gas -solid phase reaction 


] + Vw h. - W " (h - e ) 

) i—d 11 e g c 


n 

Convection 


ni 

Diffusion 


IV 


+ p C V* VT + Z p C V • VT 
& P~ _.i p. d. 


g 


V 


i 'll 

VI 


Heat Conduction Thermal Radiation Pressure Kinetic Energy 


= V . KVT 

vn 


V • q. 


R 


vm 


5 P 
3 t 

IX 


V • V • 
+ - W 


2 g 


X 
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Noting that the temperature of the gas equals that of the solid, differentiating the caloric 
equations of state, (1), (2) and (5) yields: 


3e 


at 


3e 
c 

at 


_p = 


'v at 
p 


ah 


1 = c aT 


= c 


3T 


v at 

c 


at p. at 

i 


V h. = C V T 
1 P i 


Substituting these relations into (30) and rearranging terms gives 


f 


pC +p C ■+ 
p v c v 
P c 




+ v*p. fv + v, i) + 2^P- c (v + t\ ) 

‘V ‘/J ‘ 1 p i\ V 


. VT = V ♦ KVT-7-q (31) 


+ + V • (V • T) - 


at 




The second and eighth terms can be simplified by use of the continuity equations. Using 
(14), (16) and (17), the second term becomes 


D * 0 


e W + e W + ?. h. W. 
p p c c ~ i i 
i 


Using (23), this last term is: 


Ew.i..-Z 

s 1 1 I 


W. h. + 


W h 
g g 


From (19), (20), and (21) we have: 

o • • ft 

W = -f w + w 

c c p c 

W = - (1 - f ) w + w" 

g c 7 p g 

• tt •!! 

w = - w 

g C 
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V . (V . T ) - p 


g 


Viscous Stresses 


5 

at 




XI 


( 34 ) 


Energy Equation 


Summary of Equations 


( 


pC + D C 
p v c 
P 


+ p C | 

V c g P gJ 


w 

at P 


(1 - f ) h + f e -e 
c g c c p 


\ ' . t ■ m V' _ \ 

+ / j W. h. - W (h - e ) + p C V*VT+/.p. V J h . VT 
j ii c g c' H g p ^-^i d. i 

D 11 


-i. 9 P V . V • 

= V • KVT-V* q +-^— + W + V* (V • T) 

R o t 2 g 


( 35 ) 


P g 1 ^ + 


iMV)* *•■(¥)! 


Species Continuity 


/ 8K i - 

, + V . VK( 


+7> Pi v d. 


* f 

= W. 

l 


( 36 ) 


Continuity 


ap 


g 


at + v ^ g v 


(1 - f ) w - w 

c' p c 


( 37 ) 


at p 


p = w 


( 38 ) 


A -15 


( 39 ) 



- f W + W ' 
c p c 


The various terms of this equation are identified as: 

I energy storage 

II energy absorbed due to the decomposition of the solid 

HI energy absorbed due to gas phase reactions (i. e. cracking) 

IV energy absorbed due to gas-solid phase reactions 

V energy transfer due to convection 

VI energy transfer due to diffusion 

VII energy transfer due to heat conduction 

VIII energy transfer due to thermal radiation 

IX rate of work associated with the pressure 

X kinetic energy associated with gas formation 

XI rate of work associated with the viscous stresses and kinetic energy 
The "heat of decomposition" appears in term n. 


h = (1 - f ) h + f e -e 
gf c g c c v 

P 

If the usual momentum equation could be used to simplify term XI, it would reduce to the 
familiar work of pressure forces plus the work of viscous forces (i.e. V • Vp + $ where 3> 
is the dissipation function). 



Momentum 

V 



V P 


( 41 ) 
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Diffusion (binary mixture approximation) 


p. V, = - p D 10 v K. 
i M g 12 i 


( 42 ) 


Neglecting T and q r , there are 6 + 2 (N-l) equations for the following physical variables 
(N is the total number of gaseous species) : 


T, p , p , p , P, V, V , K. 

peg d i 1 


These equations require that the following material properties (10 + 4N in all) be known 
functions of the variables. 


C , e , W 
v p P p 


C , e , f , W 
V c c c 
c 


C , h., w! , M. 

p. l l l 


K ’ T’ °12 


Discussion 


The equations developed so far represent a quite general physical model of charring 
ablation. They account for the simultaneous transfer of energy and mass within a solid 
material of variable porosity which is decomposing. The ablation gases may be flowing, 
diffusing, reacting with themselves or reacting with the char and they are not necessarily 
in local thermochemical equilibrium. 

It is generally desirable to invoke further physical assumptions in order to simplify the 
mathematical analysis and to reduce the number of required material properties, which are 
often not known. Several of these assumptions will now be discussed. 

Two approaches will be described for the simplification of the general equations derived 
above. One approach to the problem is to simplify the gas chemistry while retaining the 
gas dynamical features. The ultimate end in this approach is to assume that the gas contains 
only a single specie. Note that this assumption does not exclude gas-solid phase reactions. 
Thus, we have = 0, W ’ = 0 and the species continuity equation is superfluous. 


A-17 


Neglecting the radiant flux and using the definitions of W (19) and W (20), equations 

o 

(35) - (39) simplify to: 


Energy 


( 


p C + e C + p C 

p vp C VC g p 


g 


)*W* 


W e + W e + W, 
p p c c 


^ 

a v • KV T + jjP v W + v • (V* r) -e P- 

er ' ~ got 


at 2 g 


h h 6 ) +<p 


C V) • V T 
g p g 

(43) 


Continuity 


dP V 

rrf- + V. e V = W 

at g g 


a p 

p 

at 


w 

p 


at 


w 

c 


W + W + W = 0 
P C g 


(44) 


(45) 


(46) 


(47) 


State 



Momentum 


k 

V = - V— P 

P 


(48) 


(49) 
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The gas density used in the derivation of the equations described here is the weight of 
the gas in a given solid-gas volume divided by that volume. The gas density referred to 
within the REKAP program is the above density divided by the porosity of the material. 
Porosity is defined as: 


_ Actual Gas Volume _ Void Volume 
Total Volume ” Total Volume 

The porosity is calculated at each time step by: 



The thermal conductivity of the gas (k g ) and thermal conductivity of the solid (K ) are each 
based on their respective areas. ® 

P c is the final density of the char based on the volume of the char (i. e. if the char is 

carbon then p c is equal to the density of carbon) and p is the initial density of the virgin 

vp 

plastic before heating or any charring has taken place. 

Now simplifying the equations it will be assumed that the gas is in local thermochemical 
equilibrium and that diffusion is negligible. This means that the gas composition is a known 
function of temperature and pressure. Finally, it is assumed that radiation can be accounted 
for as an increase in effective conductivity and that the mechanical work terms are negligible 
compared to the thermal terms in the energy equation. All of these are reasonably plausible 
engineering assumptions. 

In the absence of diffusion, the species continuity equation (22) is: 


= W. 
1 


g 


SK. 
1 

at 


v • VK. 


Since the composition is a known function of temperature and pressure, this implies: 


• » 


W. = 
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+ V • VT 


+ P 
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Substituting this into the energy equation (34) and neglecting diffusion, radiation, and 
mechanical work terms we get: 



C ■ + 


v 

P 


C 


v 

c 


+ p 


gGp g) 


3T 

at 


W (1 - f ) h + f e 
p c g c c 


e 

P 



- W 


(h 


g 



5T 

a t 


v 


VT 


-i. 


g 




V T 



= V • K V T 



+ V • V P 


( 50 ) 


The continuity equations (34) - (39), thermal equation of state (40) and the momentum 
equation (41) remain unchanged. 


It is possible to simplify (50) even further by neglecting the gas density in comparison 
with the solid density, while retaining the mass flow rate term. This implies that as 

p -» 0, p V remains finite so that V -* », i. e. , the "residence time" is negligible. This 
g g 

means that the equation of state and the momentum equation are superfluous. The pressure 
is assumed to be uniform at its ambient value, which is not necessarily steady. The 
continuity equations remain unchanged except the time derivative of the gas density in (37) 
is dropped. The energy equation (50) becomes: 


( Pp 


C + p C 
v c V 
P C 


)-lf. 


. VT = V* KV T 


+ W (l-f)h + f e -e 
" 1 e g c c p 


(51) 


+ w" (h -e.) 

e g i 


This equation in one dimension is the equation solved in the REKAP program. Without 
the pressure option, and neglecting the gas-solid phase reactions and by combining equations 
(37) - (39) and (51), the following equations result: 
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Derivation 


BEKAP Program 


h + e 

* 1 - f c 0 




(1 - f ) W 
c p 


P-P, 


vp V P 


Summary of the Equations 


REKAP Equations : 
Energy: 


“p St 


( K irf “ g ( 


C +H ) f-^ + H ^ 
p g C g / Sx 


Continuity: 


backface 


M = 
g 


f _S_0_ 

J at 


Density: 


so. Ze' E/RT 

0t H °vp / 


These are the equations that are solved in the non-pressure option of the 


program. 
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Energy: 



Continuity: 

X 

P„ V = - / -^7 (p + p ) dX 

g J o t p c 

backface 

Density: 


a 

~TT (P + P ) = - a-f ) W 

at c c 7 p 


Derivation 


REKAP Program 
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+ P c 
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P 


P 


C + p 
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Pressure Option Equations: 
Energy (Gas): 


BT 


PC 


g _ 


g p at 


-pc v||. + 

g p ax 


ap a / a 

at ax K g "ax/ Q 


transferred from gas 


+ W h +-^L 
g g ax 


"if- 


p v 2 ^ 

p g ax 


w v 
g 


Energy (Solid): 


aT 


(1 - € ) p C + p C 


s a 
+ 


p r vp c v c / at ax 


aT 


(1 - € ) K 


s ax 


-( e c W 


J-Q 


W + e W / - Q . 
c vp vp J transferred from gas 


If the temperature of the gas and the temperature of the adjacent solid material are the 
same, the above two equations may be added together. 


°gV (1 - e > 


(p C + p C I 

I p V V V / 

\ vp c/J 


aT aT 
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at g g c P ' v g P g ax ax( e i. v y 
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Continuity: 



St 


w 

s 


Momentum: 


bx\ n <r P g 


which was obtained by substituting the Equation of State in the momentum equation (12). 
Boundary Conditions: 

The boundary conditions that are of concern here are those describing the material heat 
input or removal from the front and back face and the surface recession at the heated face. 
The heating of the material can be described by three methods: front face temperature (T w ), 
front face heat flux (q c and/or qw.)» an d front face convective film coefficient (q /Ah or 
VAT). c 


= W 


= W + W 
c vp 


a / g R T R 


\ 
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Each of these quantities can be a function of time. The convective film coefficient option 
is the one most commonly used for the analysis of rocket engines, however, for some 
propellant combinations, it is necessary to account for the radiation (qhgr) from the exhaust 
gases. The program includes the capability of combining the radiative flux and the convective 
heat transfer by taking a thermal balance at the front face. The thermal balance is described 
by: 





K 

w 


5 T 

ax 


The convective heat flux (q,,) is determined either from program input which is a function of 
time or it is calculated from: 




or 



q„ = 


AT 


(T - T ) 
r w 


where h^ is the recovery enthalpy and T^ is the recovery temperature. 

If the convective film coefficient is in terms of temperature rather than enthalpy, the 
specific heat (C ) of the boundary layer gases must be set equal to 1. 0 for all values of 
p bl 

gas temperature. The convective film coefficient is an input to the program and is considered 
to be a function of time. 

The fourth term in the heat balance equation is the rate of energy loss from the front 
face due to thermal radiation. It is expressed by: 


where a is the Stefan- Boltzmann constant which equals 0.476 x 10~ 12 BTU/sec Ft 2 °R 4 and e 
is the product of the surface emissivity and the configuration factor (F a ) between the point 
radiating and the cold (relative to the hot wall) external environment. 

The fifth term is the decrease in the convective heat flux due to the injection of ablation 
gases into the boundary layer. This is commonly referred to as the "blocking action effect. " 
The expressions describing the blocking action were derived from the correlation of 


A-25 



experimental data (References 14 to 34). The expressions for the blockage of the convective 
energy are: 

Laminar: 


% H 




1/3 


<f>. 


1/3 


Turbulent: 


q b = q c L 


1 - e 


• 38 <p 


] 


<p = M 

w • /A , 
q /A h 
c 


where: 

is the molecular weight of the injection gases 

Mg is the molecular weight of the boundary gases 

C T is the ratio of the specific heat of the injection gases to the specific heat of 

the boundary layer gases. C /C 

P 1 P 2 

P r is the Prandtl number of the boundary layer gases. 

• ' 

o 

M is the mass injection rate at the front face. ~ lb/sec ft 

The quantities (M 2 /M ;l ), and P r are input constants while (q c /Ah) is the value of the 
convective film coefficient. If the blocking action is expected to be significant (however, for 
most materials exposed to a rocket.engine environment, the blockage effects amount to only 
a few per-cent of the convestive heat flux) it is necessary to use the f ilm coefficient defined 
in terms of the enthalpy difference since that is how the above blocking action equations are 
correlated. Included for completeness is the laminar blocking action equation although for 
most rocket nozzle applications, the boundary layer is assumed to be turbulent. 
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The last term in the front face heat balance equation is the rate of thermal energy which 
is transferred by conduction into the material. 

The heat transfer from the back face of the material is controlled by specifying the back 
face temperature (T BF ) or heat transfer rate (q BF ) as a function of time. If radiation from 
the back face is desired, included is a routine to allow for an air-gap or a non-solid layer 
in the nozzle wall. Therefore, to account for radiation from the back face of a nozzle, the 
third layer from the last in the program is the actual nozzle backside, the second from last 
is the air-gap and the final layer is the nozzle surroundings for which the back face temperature 
is specified. Using the air-gap routine not only can radiation from the back face be accounted 
for but also natural convection and forced convection by the proper adjustment of constants. 

For the details of the mathematical equations, see Appendix B and E. 

Front Face Recession 

The front face recession is presently controlled by five methods: no melting, or recession, 
specified char length, graphite oxidation and sublimation, refrasil option and fixed melting 
temperature. The first method is normally used for the purpose of evaluating the temperature 
distribution within material which is known not to have a dimensional change. The options 
most commonly used are the specifiec char length, graphite sublimation and the fixed 
melting temperature. The fourth option which is referred to as the refrasil option is based 
on the work done by Munson and Spindler (Reference 12). For this option, the heat balance 
equation at the front face is: 


- K 


dT 


dX “ % + «hgr - V - % ■ WL,S 
where the surface recession rate is given by: 


S 



T 

w 


*2 


-V T 

e 


w 


The constants 0 lt & 2 and 0 3 are determined from experimental data. In their paper, 
Munson and Spindler listed the values of , j3 2 and j3o for silica phenolic as 0. 00917 ft/sec°R 2 
2. 0 and 1 x 10 b °R. 

The quantity (p C L), the surface or final char density and the latent heat of fusion or 
vaporization depending upon whether the material melts or is vaporized. Empirical and 
analytical (Reference 35) work done on the analysis of glassy materials within rocket 
nozzles has shown that the major portion of the surface loss to be by melting and not by 
vaporization. Therefore, the value of L for a phenolic refrasil material should be the heat 
of fusion for the char material which is primarily refrasil. The refrasil option has the 
disadvantage of being relatively slow (requires several times as much computation time as the 
fixed melting temperature option) since it must iterate on the rate of melt for each time step. 
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The fixed melting temperature option usually satisfies the rate of melt criterion after 
the first iteration. The net heat balance for the fixed melting temperature option is the 
same as for the refrasil option. However, the rate of melt (S m ) is given by: 

„ dT . 
dX ^net 

rp c L 

where p c is the density of the char and L is the latent heat of vaporization or melting 
depending on whether the material vaporizes or melts. The gasification factor T is the 
ratio of the char material which is either vaporized or melted to the total char that is lost. 
Some of the char may be lost by char popoff or some other mechanical means. The value of 
r must be determined experimentally. The front face is not allowed to recede until the front 
face temperature reaches the specified melting temperature. 

The specified char thickness option is as the name inplies, the char layer is allowed to 
grow until it reaches the specified value. Then the outer boundary moves at the same ratio 
as the reaction zone. The maximum allowable char thickness is determined by the material 
and the environment to which it is exposed. The thickness values are determined from 
experimental data. 

The graphite oxidation and sublimation option for the control of front face recession also 
accounts for the oxidation of most graphite materials including pyrolytic graphite and on 
an oxidation process, which is rate-controlled at low (1500°R) surface temperatures, but 
rapidly become diffusion-controlled as the surface temperature rises (see Figure A2 and 
A3). For the range of surface temperatures, approx, between 2500°R and 5000 R the rate 
of the overall mass loss is dominated by the slowest step, which is the counterdiffusion 
process in the multicomponent boundary layer. When the surface temperature is in this 
range, the oxidation rate levels off and becomes insensitive to the magnitude of surface 
temperature, simply because the mass loss is controlled by the diffusion of oxygen-bearing 
species to the surface rather than the specific reactivity of graphite. At even higher surface 
temperature (T w 500°R) the mass loss due to vaporization exceeds the diffusion controlled 
oxidation mass loss rate. This region is normally referred to as the sublimation regine. 

The results shown in Figure A3 were correlated (Refs. 36 and 37) and the resulting 
equations were: 


S 

m 


M 


t 


9 - 67 

M | r + 2.64 x 10 P e 
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11.05 x 10 


T 
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where the M 


is the mass loss within the diffusion controlled regime 



K 


+ K 2 


(h - C T ) 
r p bl w 


The quantities Kj and K 2 are input constants and for turbulent flow, their values for an 
air boundary layer are 4240 and 5. 77 respectively. The rate of front face recession is given 
by: 


S 

m 



surface 


The heat balance at the front face is given by: 


-K 


dT 

dX 




where 


q „ I 1 - S* (3. 96 x lo 8 ) P “* 67 e 


11. 05 x 10 4 
T 

w 


The local edge pressure P g is an input quantity which is a function of time and S* is a table 
lookup which is a function of the recovery enthalpy. 
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NOMENCLATURE 


A 

surface area 

ft 2 

C 

V 

specific heat at constant volume 

BTU/lbm °R 

C 

P 

specific heat at constant pressure 

BTU/lbm °R 

c" 

p g 

\ v . K. C . = average specific heat 
Li i pi 

BTU/lbm °R 

c 

p bl 

specific heat of the boundary layer gases 

BTU/lb 

c , 

ratio of the specific heat of the injection gases to the 
specific heat of the boundary layer gases 


D i. 

3 

D 12 

multicomponent diffusion coefficient 

ft 2 /sec 

binary diffusion coefficient 

ft 2 /sec 

e 

specific internal energy 

BTU/lbm 

®F 

energy of formation 

BTU/lbm 

~e 

unit base vector 


f 

c 

fraction of unreacted material that forms char 


h 

r 

recovery enthalpy 

BTU/lb 

h ff 

boundary layer gas enthalpy at wall temperature 

BTU/lb 

k 

permeability 

ft 2 

K 

thermal conductivity 

BTU/ft-sec ' 

K. 

1 

H. mas3 concentration 


1 

P 


K. 

l. 

3 

conductivity tensor 

BTU/ft-sec ' 

K 1 &K 2 

constants in mass loss equation 


M 

molecular weight 

lbm/mole 
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= average molecular weight 


lbm/mole 


M 

g 


M i 

M 2 


n 


P 

P 

P 

e 

Pr 

€ 


p c 

p or p 
vp vp 

6 

q c 

q hgr 

q 

rr 

q b 

q bf 


"M. 

1 

molecular weight of the injection gases 

molecular weight of the boundary layer gases 

mass injection rate at the front face 

total number of moles per unit volume; degradation 
reaction order 

pressure 

surface force per unit area 

boundary layer edge pressure 

Prandtl number of the boundary layer gases 

porosity 

final density of the char based on the volume of the char 

initial density of the virgin plastic 
convective heat flux 

heat flux due to hot gas radiation 

reradiative heat flux 

convective heat flux blocked due to mass injection 
heat flux to the backface 


qc or q 


qR or q. 


R 


R 


conduction heat flux vector 

radiant heat flux vector 
universal gas constant 


lb/lb mole 
lb/sec-ft 2 

3 

moles/ft 

lbf/ft 2 

lbf/ft 2 

lb/ft 2 


lbf/ft 3 

lbf/ft 3 

BTU/ft 2 -sec 

BTU/ft 2 -sec 

BTU/ft 2 -sec 

BTU/ft 2 -sec 

BTU/ft 2 -sec 

BTU/ft 2 -sec 

BTU/ft 2 -sec 
lbf ft/lbm mole °R 
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cH 


T 


T 

r 


T 

w 



T 


V 



temperature 
recovery temperature 
wall temperature 
temperature of backface 
stress tensor 
volume 

absolute velocity of ith species 
diffusion velocity of ith species 


R 

°R 

°R 

Ibf/ft 2 

ft 3 

ft/sec 


ft/sec 


V mass averaged velocity 


Wj or co. 


net rate of production of the ith gaseous species due 
to all chemical reactions 


. » . » 

W. or co . 
1 1 


net rate of production of ith species due to gas phase 
reactions 


W "or co” 
g g 


rate of production of gas due to gas- solid phase reactions 


X 


mole fraction 


ft/sec 


lbm/ft 3 

lbm/ft 3 

lbm/ft 3 


P 

ff. 

1 . 

3 

r 

M 

a 

T 


outward unit normal 
density 

component of the stress tensor 

viscous stress tensor 
gas viscosity 

Stefan- Boltzmann constant (0.476 x 10 * 2 Btu/sec-ft^ °R^) 
gasification ratio 


lbm/ft 3 

lbf/ft 2 

lbf/ft 2 
lbf- sec/ft 2 
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Figure Al. Schematic of a Degrading Plastic and Corresponding Profiles 
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SUBLIMATION 



Figure A2. Mass Transfer Regimes for Ablating Graphite 
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Figure A3. Normalized Ablation Rate of Graphite Over the Entire Range of Surface Temperature 




TABLE 6. GRAPHITE CLOTH/PHENOLIC RESIN - SPECIFIC HEAT 
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TABLE 7. SILICA CLOTH/PHENOLIC RESIN - THERMAL CONDUCTIVITY 
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TABLE 8. SILICA CLOTH/PHENOLIC RESIN - SPECIFIC HEAT 
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TABLE 9. PROPERTY VARIATION FOR SIGNIFICANT PROPERTIES FOR SILICA CLOTH/PHENOLIC RESIN 
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TABLE 12. NOMINAL PROPERTY VALUES FOR REKAP ANALYSIS 



SILICA CLOTH/ 
PHENOLIC RESIN 

GRAPHITE CLOTH/ 
PHENOLIC RESIN 

GRAPHITE CLOTH/ 
EPOXY RESIN 

1. Heat of Gasification 

550 

550 

550 

(Btu/lb) 

4 

4 

4 

2. Collision Frequency 

3 X 10 

3 X 10 

3 x 10 

(l/sec) 

3. Activation Energy 

Table 9 

Table 10 

Table 11 

(Btu/lb) 

4. Melting Temperature 

Table 9 

Table 10 

Table 11 

of Fibers (°R) 

5. Heat of Vapor of 

71 

Table 10 

Table 11 

Reinforced Fibers 
(Btu/lb) 

6. Wall Emissivity 

0.65 

0.8 

0.8 

7. Recovery 

Temperature (°R) 

a) N 2 O 4 /Aerozine 50 

4565 



b) of 2 /b 2 h 6 

— 

6480 

6480 

8. Film Coefficient 
(Btu/ft2sec°R) 

a) 1. 2 in. dia. Throat 

0.294 

0.425 

0.425 

b) 7. 82 in. dia. Throat 

0.223 


— 

9. Specific Heat of 

0.75 

0.75 

0.75 

Ablation Gases 
(Btu/lb°R) 

10. Molecular Weight of 

30 

30 

30 

Ablation Gases 
ll.Virgin Plastic Density 

Table 9 

Table 10 

Table 11 

(lb/ft 3 ) 

12. Char Density (lb/ft 3 ) 

Table 9 

Table 10 

Table 11 

13.Thermal Conductivity 

Figure 4 

Table 10 

Table 11 

(Btu/ft. sec°F) 
14.Specific Heat 

(Nom. Curve) 
Table 9 

Figure 85 

Figure 85 

(Btu/lb) 


(Nom. Curve) 

(Nom. Curve) 

15,Order of Reaction 

2 

2 

2 
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TABLE 11. PROPERTY VARIATION FOR SIGNIFICANT PROPERTIES FOR GRAPHITE CLOTH/EPOXY RESIN 





CONCLUSIONS AND RECOMMENDATIONS 


The surface recession of the silica phenolic material is most affected by the melting 
temperature of reinforcing fibers. The recession rate was decreased by a factor of three 
by increasing the melting temperature from 3000°F to 4000OR. Increasing the virgin plas- 
tic density from its minimum to its maximum value, increasing the activation energy from 
its minimum to its maximum value and increasing the specific heat from its minimum to 
its maximum value causes only a 91 per cent change in the surface recession rate of silica 
phenolic. Therefore, the greatest improvement in this material would be produced by in- 
creasing the fiber melting temperature. 

The surface recession rates of the graphite phenolic and graphite epoxy are most af- 
fected by the surface reaction constants (measure of reactivity of the surface material with 
the boundary layer gases). The surface recession rate of the graphite phenolic material is 
increased by a factor of 8 in going from the minimum to the maximum values of the surface 
reaction constants. For graphite epoxy, changing the surface reaction constants from their 
minimum to maximum values increases the recession rate by a factor of ten. Allowing the 
three material properties (char density, thermal conductivity, virgin material density) 
which have the next strongest influence on the surface recession rate to go to their extreme 
values causes only a 44 per cent change in the recession rate of graphite epoxy. Therefore, 
to minimize the surface recession of the graphite materials, additives or methods of pro- 
tecting the nozzle walls from the reactive chemical species in the boundary layer are re- 
quired. 
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Figure 3. Preliminary Thermal Conductivity of Graphite Cloth/Epoxy Resin for Virgin and Charred Condition 
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Figure 4. Preliminary Specific Heat of Silica Cloth/Phenolic Resin for Virgin and Charred Condition 
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for Virgin and Charred Condition 




RECOVERY TEMPERATURE, DEG R 

Figure 7. Surface Recession Rate Versus Recovery Temperature 
Silica Cloth/Phenolic Resin 
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Figure 8. Surface Recession Rate Versus Film Coefficient Silica Cloth/Phenolic Resin 
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Figure 9. Surface Recession Rate Versus Virgin Density Silica Cloth/Phenolic Resin 
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Figure 16. Surface Recession Rate Versus Specific Heat of Ablation Gases 

Silica Cloth/Phenolic Resin 
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Figure 19. Surface Recession Rate Versus Recovery Temperature 
Graphite Cloth/Phenolic Resin 
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Figure 20. Surface Recession Rate Versus Surface Reaction Rate 
Constant K Graphite Cloth/Phenolic Resin 
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Figure 21 . Surface Recession Rate Versus Film Coefficient 
Graphite Cloth/Phenolic Resin 
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Figure 23. Surface Recession Rate Versus Char Density 
Graphite Cloth/Phenolic Resin 
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Figure 24. Surface Recession Rate Versus Thermal Conductivity Ratio 
Graphite Cloth/Phenolic Resin 
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Figure 25. Surface Recession Rate Versus Virgin Density 
Graphite Cloth/Phenolic Resin 
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Figure 27. Surface Recession Rate Versus Specific Heat Ratio 
Graphite Cloth/Phenolic Resin 
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Figure 28. Surface Recession Rate Versus Activation Energy 
Graphite Cloth/Phenolic Resin 
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Figure 31. Surface Recession Rate Versus Heat of Gasification 
Graphite Cloth/Phenolic Resin 
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Figure 33. Surface Recession Rate Versus Surface Reaction Rate 
Constant K-^ Graphite Cloth/Epoxy Resin 
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Figure 38. Surface Recess 
Graphite Clc 
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Figure 39. Surface Recession Rate Versus Activation Energy 
Graphite Cloth/Epoxy Resin 
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Figure 40. Surface Recession Rate Versus Specific Heat Ratio 
Graphite Cloth/Epoxy Resin 
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Figure 42. Surface Recession Rate Versus Collision Frequency 
Graphite Cloth/Epoxy Resin 
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Figure 43. Surface Recession Rate Versus Wall Emissivity 
Graphite Cloth/Epoxy Resin 
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Figure 45. Mass Loss Rate Versus Char Density Silica Cloth/Phenolic Resin 
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Figure 46. Mass Loss Rate Versus Melting Temperature 
of Reinforcing Fibers 
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Figure 47. Mass Loss Rate Versus Activation Energy Silica Cloth/ Phenolic Resin 
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Figure 51. Mass Loss Rate Versus Thermal Conductivity Ratio 
Silica Cloth/ Phenolic Resin 
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Mass Loss Rate Versus Heat of Vaporization 
Silica Cloth/Phenolic Resin 
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Figure 57. Mass Loss Rate Versus Char Density 
Graphite Cloth/Phenolic Resin 
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Figure 61. Mass Loss Rate Versus Thermal Conductivity Ratio 
Graphite Cloth/ Phenolic Resin 
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Figure 64. Mass Loss Rate Versus Surface Reaction Rail Constant K 

Graphite Cloth/Phenolic Resin 
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Figure 65. Mass Loss Rate Versus Specific Heat Ratio 
Graphite Cloth/ Phenolic Resin 
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Figure 66. Mass Loss Rate Versus Specific Heat of Ablation Gases 
Graphite Cloth/Phenolic Resin 
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Figure 67. Mass Loss Rate Versus Heat of Gasification 
Graphite Cloth/Phenolic Resin 
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Figure 68. Mass Loss Rate Versus Wall Emissivity 
Graphite Cloth/Phenolic Resin 
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Figure 69. Mass Loss Rate Versus Char Density 
Graphite Cloth/Epoxy Resin 
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Figure 71. 
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Figure 72. Mass Loss Rate Versus Thermal Conductivity Ratio Graphite Cloth/Epoxy Resin 
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Mass Loss Rate Versus Activation Energy Graphite Cloth/Epoxy Resin 
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Figure 78. Mass Loss Rate Versus Specific Heat of Ablation Gases 
Graphite Cloth/Epoxy Resin 
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Figure 82. Variations in Thermal Conductivity Due to Changes in Lamination Angle, Resin Content, and Lot-to-Lot 

Variation of Graphite Cloth/Phenolic Resin (Literature Survey) 
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Figure 84. Thermal Conductivity of Graphite Cloth/ Phenolic Resin 

(Literature Survey) 
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Figure 85. Specific Heat of Graphite Cloth/Phenolic Resin and Graphite Cloth/Epoxy Resin 

(Literature Survey) 





Figure 86. Variations in Thermal Conductivity Due to Changes in Lamination Angle, Resin Content, and Lot -to- Lot 

Variation of Silica Cloth/Phenolic Resin (Literature Survey) 
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of Silica Cloth/Phenolic Resin (Literature Survey) 



NOTE: THE DASHED LINES BOUND THE REGION ±50% OF THE SOLID LINE 
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Figure 88. Thermal Conductivity of Silica Cloth/Phenolic Resin (Literature Survey) 



NOTE: THE DASHED LINES BOUND THE REGION ±15% OF THE SOLID LINE 
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Figure 89. Specific Heat of Silica Cloth/Phenolic Resin (Literature Survey) 
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Figure 90. Thermal Conductivity of Graphite Cloth/Epoxy Resin (Literature Survey) 
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Figure 91. Surface Recession Rate Versus Virgin Plastic Density Silica 
Cloth/Phenolic Resin A (1 . 2 Inch Diameter Throat) 


123 




SURFACE RECESSION RATE ~ IN. /SEC 


0.040 


0.036 


0.032 


0.028 


0.024 


0.020 


0.016 


0.012 


0.008 


0.004 


0 


MELTING TEMPERATURE 


3000°R 



= 3500 R 

ACTIVATION ENERGY = 75, 600 BTU/LB 
|” ACTIVATION ENERGY = 48, 600 BTU/LB 
ACTIVATION ENERGY = 21, 600 BTU/LB- 



ACTIVATION ENERGY = 75, 600 BTU/LB 
ACTIVATION ENERGY = 48, 600 BTU/LB 
ACTIVATION ENERGY = 21, 600 BTU/LB 

MELTING TEMPERATURE 


|— GENERAL NOTE: 

SPECIFIC HEAT RATIO 

WHEN VIRGIN PLASTIC DENSITY =88 LB/FT U ,CHAR DENSITY= 85 LB/FT C 


= 103 LB/FT' 5 
= 118 LB/FT 3 


= 92 LB/FT 3 
= 99 LB/FT 3 


82 88 94 100 106 112 

VIRGIN PLASTIC DENSITY ~ LB/FT 3 


118 


124 


Figure 92. Surface Recession Rate Versus Virgin Plastic Density Silica 
Cloth/Phenolic Resin A (1 . 2 Inch Diameter Throat) 
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Figure 93. Surface Recession Rate Versus Virgin Plastic Density Silica 
Cloth/Phenolic Resin A (1.2 Inch Diameter Throat) 
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Figure 94. Surface Recession Rate Versus Specific Heat Ratio Silica Cloth/ 
Phenolic Resin A (1. 2 Inch Diameter Throat) 
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Figure 95. Surface Recession Rate Versus Specific Heat Ratio Silica Cloth/ 
Phenolic Resin A (1. 2 Inch Diameter Throat) 
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Figure 96. Surface Recession Rate Versus Specific Heat Ratio Silica Cloth/ 
Phenolic Resin A (1. 2 Inch Diameter Throat) 
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Figure 97. Surface Recession Rate Versus Activation Energy Silica 
Cloth/Phenolic Resin A (1. 2 Inch Diameter Throat) 
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Figure 98. Surface Recession Rate Versus Activation Energy Silica 
Cloth/ Phenolic Resin A (1.2 Inch Diameter Throat) 
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Figure 99. Surface Recession Rate Versus Activation Energy Silica 
Cloth/Phenolic Resin A (1. 2 Inch Diameter Throat) 
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Figure 100. Surface Recession Rate Versus Virgin Plastic Density Silica Cloth/ 
Phenolic Resin B (7. 82 Inch Diameter Throat) 
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Figure 101. Surface Recession Rate Versus Virgin Plastic Density Silica Cloth/ 
Phenolic Resin B (7. 82 Inch Diameter Throat) 
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Figure 102. Surface Recession Rate Versus Virgin Plastic Density Silica Cloth/ 
Phenolic Resin B (7. 82 Inch Diameter Throat) 
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Figure 103. Surface Recession Rate Versus Specific Heat Ratio 
Silica Cloth/Phenolic Resin B 
(7. 82 Inch Diameter Throat) 
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Figure 104, Surface Recession Rate Versus Specific Heat Ratio 
Silica Cloth/Phenolic Resin B 
(7. 82 Inch Diameter Throat) 
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Figure 105. Surface Recession Rate Versus Specific Heat Ratio 
Silica Cloth/Phenolic Resin B 
(7. 82 Inch Diameter Throat) 
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Figure 106. Surface Recession Rate Versus Activation Energy 
Silica Cloth/Phenolic Resin B 
(7. 82 Inch Diameter Throat) 
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Figure 108. Surface Recession Rate Versus Activation Energy 
Silica Cloth/Phenolic Resin B 
(7. 82 Inch Diameter Throat) 
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Figure 107. Surface Recession Rate Versus Activation Energy 
Silica Cloth/Phenolic Resin B 
(7. 82 Inch Diameter Throat) 
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Figure 109. Surface Recession Rate Versus Char Density 
Graphite Cloth/Phenolic Resin 
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Figure 110. Surface Recession Rate Versus Char Density 
Graphite Cloth/Phenolic Resin 
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Figure 111. Surface Recession Rate Versus Char Density 
Graphite Cloth/Phenolic Resin 
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Figure 112. Surface Recession Rate Versus Virgin Plastic Density 
Graphite Cloth/Phenolic Resin 
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Figure 113. Surface Recession Rate Versus Virgin Plastic Density 
Graphite Cloth/Phenolic Resin 
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Figure 114. Surface Recession Rate Versus Virgin Plastic Density 
Graphite Cloth/Phenolic Resin 
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Figure 115, Surface Recession Rate Versus Thermal Conductivity Ratio 

Graphite Cloth/Phenolic Resin 
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Figure 116. Surface Recession Rate Versus Thermal Conductivity Ratio 

Graphite Cloth/Phenolic Resin 
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Figure 118. Surface Recession Rate Versus Char Density 
Graphite Cloth/Epoxy Resin 


150 



SURFACE RECESSION RATE- IN. /SEC 



CHAR DENSITY - LB/FT 3 


Figure 119. Surface Recession Rate Versus Char Density 
Graphite Cloth/Epoxy Resin 
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Figure 120. Surface Recession Rate Versus Char Density 
Graphite Cloth/Epoxy Resin 
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Figure 121. Surface Recession Rate Versus Virgin Plastic Density 
Graphite Cloth/Epoxy Resin 
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Figure 122. Surface Recession Rate Versus Virgin Plastic Density 
Graphite Cloth/Epoxy Resin 
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Figure 123. Surface Recession Rate Versus Virgin Plastic Density 
Graphite Cloth/Epoxy Resin 
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Figure 124. Surface Recession Rate Versus Thermal Conductivity Ratio 

Graphite Cloth/Epoxy Resin 
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Figure 125. Surface Recession Rate Versus Thermal Conductivity Ratio 

Graphite Cloth/Epoxy Resin 
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Figure 126. Surface Recession Rate Versus Thermal Conductivity Ratio 

Graphite Cloth/Epoxy Resin 
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INTRODUCTION 


In this appendix, a thermal ablation model is derived for a thermosetting plastic. Con- 
sideration is given first to the general three-dimensional case. Simplifications are then 
introduced to obtain an equation which reasonably satisfies the physical model. 

The philosophy of this derivation is to start from fundamental physical principles and to 
utilize the concepts of continuum mechanics to proceed in a step-by-step fashion, listing 
all assumptions. 

Figure A1 shows a cross-section of the ablation model. Initially, the outer boundary 
coincides with the broken line as indicated. The ambient temperature is low enough so that 
no chemical reactions occur within the plastic. Furthermore, the outer boundary tempera- 
ture is the same as its surroundings and, therefore, radiation to or from the front face 
is zero. 

Convective and radiative heat fluxes (arbitrary with time) are impressed on the outer 
boundary. As a consequence of thermal conduction, laminates of the plastic near the outer 
boundary increase in temperature and the front face begins to radiate heat. In time, the 
hotter laminates undergo a chemical reaction which converts the virgin plastic into hydro- 
carbon gas and a porous char residue. 

The gas pressure within the porous char increases as the virgin material undergoes 
chemical reaction. As a consequence, a pressure profile is established throughout the 
porous region causing the gas to flow to adjacent pores of lower pressure. In general, the 
gas flow will be to the outer boundary and result in thermal energy being introduced due to 
friction. Heat transfer will occur between char and gas if their respective temperatures are 
different. Varying temperature or pressure changes, or any combination of these two con- 
ditions, can result in chemical changes in the gas (cracking or recombination), which will 
absorb or generate thermal energy. As the gas passes the outer boundary, a portion of the 
convective heat flux is blocked. As more and more heat enters the front face, reacting 
laminates will completely de-gas, thus forming a char layer while moving the reaction zone 
deeper into the body. And, of course, the outer boundary moves as a result of structural 
failure, oxidation, or both. If the outer boundary temperature becomes high enough, the 
char layer will either melt as in the case of the material having silica fibers, or undergo 
surface reaction with the boundary layer gases as for the graphite materials. 

Physical Model 

The physical model is that of a multicomponent flow of chemically reacting gases through 
a porous media which is itself undergoing chemical reactions. The ablation material consists 
of unreacted solid (denoted by subscript p), which decomposes to a porous solid (subscript C) 
and gaseous products of reaction (subscript g). The decomposition process can be schemati- 
cally represented as: 

P(S) - C(S) + G(g) 
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Before decomposition begins, the ablation material consists solely of unreacted solid. After 
the process has gone to completion, only solid and gaseous products of reaction exist. 

All densities are based on the same unit reference volume of the mixture (solid and gas). 
Consequently, as the decomposition proceeds at a given location, p i decreases from some 
initial value to zero while p c is simultaneously increasing from zero to some final value. 

The gaseous ablation products are formed by the decomposition of the unreacted solid 
material. They are a mixture of many different chemical species which flow and diffuse 
through the porous solid. The various species may react with one another in the gas phase 
resulting in the familiar ’’cracking" effect. They may also react with the surrounding solid 
material, causing a reduction (or increase) in solid density. 

In order to validly apply continuum theory to a porous media, all quantities are pre- 
sumed to be suitably averaged over a small area - small with respect to the macroscopic 
dimensions of the material but large with respect to pore size. It is assumed that the ratio 
of pore area to total area is the same as that of pore volume to total volume, the latter 
quantity being the definition of porosity. 

The solid species remain stationary as the displacements due to thermal expansion, a 
stress field and/or changes in molecular structure are generally negligible. All species 
are considered to be pure substances. External body forces (e.g. gravity) have been neg- 
lected as they are small for all practical applications. 


Equations of State 

The caloric equation of state for each solid specie is assumed to be of the form: 


e = e (T) 
P P 


Thus, for any process: 


e =e (T) 
c c 



( 1 ) 


( 2 ) 


The internal energy accounts for thermal and chemical energy. The solid species do not have 
a thermal equation of state as their densities are determined by the application of non- 
equilibrium reaction kinetics. 
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The gaseous products are assumed to be a mixture of chemical reacting perfect gases. 
Thus, the thermal and caloric equations of state for each specie are: 



Note that P. and p i are partial quantities which are based on a reference volume of the 
entire mixture (solid plus gas). 

For the gaseous mixture as a whole, we have (assuming Dalton's Law of Partial Pres- 
sures is valid): 


P = p 



M = 
g 



e 

g 


h 

g 


■Evi 


■2>t h i 



( 6 ) 


( 7 ) 


( 8 ) 

0 ) 


Note that these assumptions imply that pressure, stress, chemical reactions, etc. have a 
negligible effect on the specific internal energy of each species. Obviously, they do affect 
the amount of each species present at a given location and thus they do effect the total 
energy. 



Diffusion Velocities 


In the flow of multicomponent gases, diffusion currents are generated by gradients in 
concentration, pressure and temperature. For the present problem, pressure and thermal 
diffusion effects should be small and so they are neglected. The velocity of the ith species 
relative to a fixed coordinate system is defined as V.. The mass-averaged or observable 
velocity of the total gas flow is defined as: 


V = 


g 


E 

i 


P- V. 
1 1 


The diffusional velocity of the ith species (V d .) is defined as the velocity of the ith specie 
relative to the mass-averaged velocity. 



Note that: 


V? 'iV? » i * 





p v = o 

g 


To summarize, the absolute velocity of the ith species is given by the vector sum of the 
mean flow velocity and the diffusional velocity of the ith species, and the mass-averaged 
diffusional velocity is zero. 


V. = V + V ’ 
i d. 

l 

_ ( 10 ) 

£ pv d =0 

i i 


For ^ordinary concentration diffusion in a multicomponent gas, a first order approximation 
for V d . is that it depends linearly upon the concentration gradients of all species. For a 
mixture of perfect gases (p. 569, Reference 4): 


p. V . = 
i d. 

l 


g 


y. M. M. P. v X. 

tt j 1 1 1 
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Use of this equation results in a formidable mathematical problem to determine the compo- 
sition of the mixture. Also, since we are dealing with transport phenomena in a porous 
media, its accuracy is not assured. It has been noted by Von Karman (Ref. 1) that ". . . 
the process in a multicomponent mixture is so complicated that one mostly uses an approxi- 
mation by considering the diffusion between one appropriately chosen component and the 
mixture of the rest replaced by a homogeneous gas of average characteristics, " i.e. , an 
effective binary mixture insofar as diffusion is concerned. 

With this approximation, the diffusion velocity is related to the mass concentration by 
Fick's Law: 

P iV - p g D 12 VK i (n > 

The concept of an effective binary mixture would be a useful starting point in accounting for 
the effects of diffusion. Probably the largest error in this approximation is that the diffusion 
coefficient for each specie is the same. 

Momentum Equation 

Experimental evidence for the flow of a gas through a porous media indicates that the 
usual momentum equation of fluid mechanics does not apply (e.g. Reference 9). Conse- 
quently, it must be replaced by an empirical relationship between velocity and pressure. 

For the flow of a homogeneous gas through a porous media at low velocities, Darcy's Law 
is reasonably accurate (Reference 9). Very little is known about the present case of chemi- 
cally reacting flow through a media of variable porosity. It will be assumed that Darcy's 
Law gives an adequate representation for the present problem, although other forms could 
be used if desired. Thus: 



where k is the permeability of the charring material and ^ is the viscosity of the ablation 
gases. These quantities are normally determined by experiments. 

Continuity Equations 

The principal of conservation of mass as applied to the ith gaseous specie within a 
stationary control volume says that the rate at which mass is accumulated within the volume 
equals the rate at which mass is transported out by convection and diffusion plus the net rate 
of production due to chemical reaction. The mathematical statement of this is: 



( 13 ) 
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